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Abstract 

Taking the two-dimensional Ising model for example, short-time 
behavior of critical dynamics with a conserved order parameter is in- 
vestigated by Monte Carlo simulations. Scaling behavior is observed, 
but the dynamic exponent z is updating schemes dependent. 
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Traditionally, it is known that there exists universal dynamic scaling be- 
havior in the long-time regime of critical dynamic evolution, since both the 
spatial correlation length and the correlation time are very large. In recent 
years, however, it is discovered that universal scaling behavior emerges al- 
ready in the macroscopic short-time regime, after a microscopic time scale 
tmic which is large enough in the microscopic sense ^, |^, H, ||, |], |^. Im- 
portant is that one should take care of the macroscopic initial conditions 
carefully. 

For a random initial state, for example, at the critical point the non- 
equilibrium spatial correlation length increases in time by a power law t^^^ 
due to the divergent correlation time 0]. This implies a kind of memory 
effect in the system. Even though the spatial correlation length (~ t^/^) is 
still small in the short-time regime, one may expect dynamic scaling. For a 
small initial magnetization, an extra critical exponent is sufficient to describe 
its scaling behavior ^, 0]. In general, one needs a characteristic function 
for an arbitrary magnetization P, ^. If the initial spatial correlation length 
is non-zero, it induces corrections to scaling |^. 

Since 1989, analytical calculations for short-time dynamic scaling behav- 
ior based on renormalization group methods 11] have extended to dif- 
ferent dynamics classified by Hohenberg and Halperin [Q. However, Monte 
Carlo simulations are still limited in dynamics of model A [0, |1^, |15|, |16 



17| , p!8| , p!9| , |2^, |2T|]. Actually, due to severe critical slowing down, Monte Carlo 



simulations for long-time behavior of critical dynamics beyond model A are 
also in preliminary stage. On the other hand, short-time dynamic scaling is 
not only conceptually interesting but also practically important. It provides 
new methods for numerical measurements of not only dynamic exponents 
but also static exponents as well as the critical temperature [p2| , These 



methods do not suffer from critical slowing down. Therefore it is important 
to investigate short-time scaling behavior of critical dynamics beyond model 
A with Monte Carlo methods. 

Dynamics of model A is a kind of relaxation dynamics without relevant 
conserved quantities |T^. If we consider only dynamic relaxation processes. 



dynamics with different relevant conserved quantities is classified into model 
B, C and D. In this letter, taking the two-dimensional Ising model for exam- 
ple, with Monte Carlo methods we study dynamics with a conserved order 



parameter, which is called model B, according to Ref. [|12 



In dynamics of model B, to guarantee criticality the conserved order pa- 
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rameter must be set to zero. The equilibrium state of model B is in a same 
universality class of model A. For lattice models, the critical temperature in 
dynamics of model B is also the same as that of model A. Of course, the 
dynamic exponent z of model B is different from that of model A. For the 0^ 
theory with a generalized Langevin equation in which the order parameter is 
conserved, in the long-time regime of dynamic evolution it has been derived 



that 2 = 4 — 1], with r] being the well known static exponent |T2|. In Ref. 
[||], it is argued that there exists also short-time dynamic scaling and no new 
independent critical exponents emerge. 

For the Ising model, a realization of dynamics of model B is considered 
to be Kawasaki dynamics in which two neighboring spins exchange their 
values in a flip. Early simulations with small lattices for dynamic behavior 
in equilibrium yield a dynamic exponent z consistent with the theoretical 



prediction z = A — rj [^|. In this paper, we study the non-equilibrium 
dynamic process starting from a random initial state. We measure the spatial 
correlation function 

C{r,t) = l-^Y.{S^it)S^+r{t)) • (l) 

Here Si is a Ising spin, d = 2 is the space dimension and L is the lattice size. 
In updating, we randomly pick up neighboring spin pairs. When L*^ pairs 
of spins are updated, we define one Monte Carlo time step. At the critical 
point, standard scaling form for C{r,t) is written as 

C{r,t)=t-'^/'^^F{^). (2) 

Here (3 and u are static exponents and z is the dynamic exponent. Usually 
such a scaling form is assumed to be valid already in the macroscopic short- 
time regime, after a time scale tmic which is large enough in the microscopic 
sense. If a Monte Carlo time step is considered as a typical microscopic time 
unit, tfnic may be several to some hundred time steps. For dynamics of model 
A, this is indeed the case 0. 

In Fig. |I|, the spatial correlation function C{r,t) for Kawasaki dynamics 
with the heat-bath algorithm is displayed. The lattice size is L = 512 and 
samples of initial configurations for averaging are 100. Maximum updating 
time is 64000. Here finite size effect is apparently negligible. For example, 
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C{r,t) at r = 80 is already very small. Extra simulations up to t = 256000 
show similar behavior as in Fig. |l], but one needs bigger lattices. 

Carefully looking at Fig. |I], we conclude that there is no scaling behavior 
as described by Eq. (0) , even though it seems there is a kind of self-similarity 
during the dynamic evolution. Let us concentrate our attention on the first 
peak in r direction (for a fixed time t). If there would be scaling, the first 
peak for different time t's should overlap each other after properly rescaling 
r and C(r, t) according to Eq. (||). We remind ourselves that the exponent 
2(3 /vz is positive. Therefore the peak for a larger t should be lower than that 
for a smaller t. But it is not the case in Fig. |1]. This fact indicates that at 
least for the short-time behavior, Kawasaki dynamics of the Ising model can 
not be described by the generalized Langevin equation of model B for the 0^ 
theory given in Ref. . 



Why there exists no short-time dynamic scaling behavior in Kawasaki 
dynamics is not very clear for us. What we observe is that in Kawasaki 
dynamics the order parameter is conserved very locally by exchanging two 
neighboring spins. Spin exchanging takes place only along domain walls. 
Smaller domains can not be so easily created during time evolution. As a 
result, the distribution of domain sizes (at a fixed time t) is rather sharp. 
This is the origin of the oscillating mode in C (r, t) as shown in Fig. |I[ This 
oscillating mode violates the scaling form. 

Dynamics with a conserved order parameter is not unique. For example, 
we may consider an alternative dynamic process in which the order parameter 
is only globally conserved. A simple case is that we exchange two randomly 
separated spins. In this case, scaling behavior is observed. To show the scal- 
ing behavior is rather robust, we have performed the simulations with dif- 
ferent updating schemes and algorithms. In Fig. ^ (a), correlation functions 
C{r,t) with the heat-bath algorithm are displayed for two different updating 
schemes. Solid lines are from a scheme in which one spin sweeps over the 
lattice regularly, while dashed lines are from a scheme in which both spins 
are randomly selected. Total samples of initial configurations for averaging 
are 620 and 220 respectively. Obviously, both updating schemes yield almost 
the same results, even though some fluctuation is seen for larger times. In 
Fig. 1^ (b), similar curves are seen for both the heat-bath and the Metropolis 
algorithms. But Metropolis dynamics evolves slightly faster. Important in 
Fig. H is that the oscillating mode is swept away. Comparing Fig. |l] and H, 
we observe that Kawasaki dynamics is much slower. 
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Now we study scaling behavior of C{r,t). According to Eq. (Q), if the 
dynamic exponent z and the static exponent (3/i> are properly chosen, C(r, t) 
at different time t's collapse on a single curve. In other words, from scaling 
collapse of C(r, t) one can estimate z and P/u. Very careful analysis shows 
that the microscopic time scale tmic's for the heat -bath and the Metropolis 
algorithms are less than 320 and 160 respectively. Also, up to the time 
t = 5120 there exists no visible finite size effect for both algorithms with 
L = 512. In Fig. ^, scaling plot for C{r,t) is shown. All data points 
collapse nicely. The fitted exponents are z = 2.315(5) and 2[3/u = 0.252(3) 
for the heat-bath algorithm and z = 2.330(15) and 2/5/z/ = 0.251(2) for the 
Metropolis algorithm. Results for both algorithms are well consistent. The 
values for the exponent 2/3/i/ agree remarkably with the exact value 0.25 
obtained in equilibrium. These results strongly support the scaling form in 
Eq. d). 

If we fix the static exponent 2(3 /v to its exact value in the scaling collapse, 
the resulting dynamic exponent zis z = 2.320(5) for the heat-bath algorithm 
and z = 2.335(10) for the Metropolis algorithm. Taking into account all 
above measurements, we conclude 

z = 2.325(10) . (3) 

This value of z is significantly different from z = 4 — 7/ = 3.75. In other 
words, the dynamics is not in a same universality class of the 0^ theory with 
the generalized Langevin equation in Ref. [|12|. The value z = 2.325(10) is 



relatively close to 2; = 2.165(10) for dynamics of model A [0, but also clearly 
different. 

Another interesting dynamic observable is the auto-correlation function 

m = Y,T.{s^{o)sm ■ (4) 

In Refs. p5|, [l|, it has been derived that at the critical point, for dynamics 
of model A, A{t) decays by a power law 

A{t) ~ . (5) 

Here the exponent A = d/ z—{2—z—ri—riQ/2)/ z, andr^o is an exponent induced 
from the divergence related to the initial condition |25, |l|]. For dynamics of 
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model B for the theory discussed in Ref. |]12|, it is pointed out in Ref. 
rjo = 0, since there is no divergence related to the initial condition. 

In Fig. ^, the auto-correlation A{t) in our simulations for Fig. ^is plotted 
in double-log scale. Power law behavior is observed starting around a time 
t = 20. After t = 1000, fluctuation becomes large. In the time interval 
[20, 1000], we measure the exponent A = 0.675(3) and 0.679(13) for the heat- 
bath and the Metropolis algorithms respectively. As an average, we estimate 

A = 0.677(8) . (6) 

Obviously, it does not yield rjo = and also is different from A = 0.737(1) 
for dynamics of model A [j^. 

In conclusions, we have investigated short-time behavior of critical dy- 
namics for the two-dimensional Ising model with a conserved order parame- 
ter. For a locally conserved order parameter (Kawasaki dynamics), standard 
scaling behavior does not exist. If we release the condition such that the or- 
der parameter is only globally conserved, universal scaling behavior emerges. 
However, the dynamics is not in a same universality class of the 0^ the- 
ory with the generalized Langevin equation which is called model B in Ref. 
[|1^. Very probably both z and A are independent critical exponents in this 
case. It is actually interesting how to construct a corresponding continuum 
model and solve it with renormalization group methods. It seems that there 
might be some similarity at certain points between this dynamics and kinetic 
roughening p6[ . 



Erratum 

We regret that the data for Fig. |^ are obtained with an incorrect updating 
scheme. The discussions and conclusions based on Fig. |I| must be modified. 
A correct version of Fig. |I| is displayed in Fig. |^. The scaling behavior 
described by Eq. (|l]) is observed. The dynamic exponent z depends slightly 
on the time t. From our data, we could only conclude that z is 3.95(10). This 
is somewhat bigger than the theoretical value z = 4 — t] = 3.75 calculated 
with the 0^ theory. Probably there exist still corrections to scaling. To 
perform simulations to larger time t's and remove these corrections to scaling. 
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however, is not so easy since z has a big value. From the auto-correlation, we 
estimate the exponent A = 0.495(10). It is consistent with the scaling relation 
\ — d/z derived with the (f)^ theory (formally A = d/ z — {A — z — r] — riQ/2)/ z — 
d/ z + riQ/2z, but t^q = 0). 
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Figure 1: Correlation functions C(r, t) for Kawasaki dynamics with the heat- 
bath algorithm. The lattice size is L = 512. From left to right, curves 
correspond to the time t = 2000, 4000, 8000, 16000, 32000 and 64000. 
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Figure 2: Correlation functions C(r, t) for dynamics with a globally conserved 
order parameter. The lattice size is L = 512. (a) The heat-bath algorithm. 
Solid lines are from a updating scheme in which one spin sweeps over the 
lattice regularly. Dashed lines are from a scheme in which both spins are 
randomly selected, (b) The solid lines are the same as in (a) but dashed 
lines are obtained with the Metropolis algorithm. In both (a) and (b), from 
left to right, curves correspond to the time t = 80, 160, 320, 640, 1280, 2560 
and 5120. 
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Figure 3: Scaling plot for C{r,t). (a) The solid line is for ti = 320, while 
o, □, A and * fitted to the curve of ti = 320 are those for t = 640, 1280, 
2560 and 5120 but r and C are rescaled according to Eq. (H) to r/^t/tiY^^ 
and (t/ti)2/^/''^C with z = 2.315 and 2p/iy = 0.252. (b) The sohd line is for 
ti = 160, while o, □, A, * and o are for t = 320, t = 640, 1280, 2560 and 
5120 with z = 2.330 and 2/5/z/ = 0.251. 
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Figure 5: Correlation functions C{r, t) for Kawasaki dynamics with the heat- 
bath algorithm. The lattice size is L = 512. From left to right, curves 
correspond to the time t = 2000, 4000, 8000, 16000, 32000 and 64000. 
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